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Abstract 

Analytical formulas for next-to-leading order temperature corrections to the 
thermal state variables of interacting nucleons in bulk matter are derived in the 
degenerate limit. The formalism developed is applicable to a wide class of non- 
relativistic and relativistic models of hot and dense matter currently used in 
nuclear physics and astrophysics (supernovae, proto-neutron stars and neutron 
star mergers) as well as in condensed matter physics. We consider the general 
case of arbitrary dimensionality of momentum space and an arbitrary degree 
of relativity (for relativistic mean-field theoretical models). For non-relativistic 
zero-range interactions, knowledge of the Landau effective mass suffices to com¬ 
pute next-to-leading order effects, but in the case of finite-range interactions, 
momentum derivatives of the Landau effective mass function up to second or¬ 
der are required. Numerical computations are performed to compare results 
from our analytical formulas with the exact results for zero- and finite-range 
potential and relativistic mean-field theoretical models. In all cases, inclusion 
of next-to-leading order temperature effects substantially extends the ranges of 
partial degeneracy for which the analytical treatment remains valid. 
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1. INTRODUCTION 

Homogeneous bulk matter comprised of fermions is commonly encountered 
in astrophysics, condensed matter physics, and nuclear physics. For extreme de¬ 
generate to near-degenerate conditions which prevail when the temperature is 
5 small compared to the Fermi temperature, Landau’s Fermi Liquid Theory (FLT) 
has been a useful guide to describe the thermodynamic and transport proper¬ 
ties of the system (see, e.g., [T] and references therein). The equation of state 
(EOS) of dense matter in cold and catalyzed neutron stars, for example, is dom¬ 
inated by the zero-temperature properties (which predominantly determine the 
10 structure of and neutrino interactions within the star) while hnite-temperature 
corrections (important for the cooling of neutron stars) are adequately given by 
the degenerate limit expressions from FLT. The leading order FLT corrections 
to the energy density and pressure are quadratic in the temperature; corrections 
to the entropy and specific heats are linear in the temperature. However, matter 
15 in supernovae and proto-neutron stars [30, especially in situations in which 
collapse to a black hole occurs, may reach temperatures exceeding the Fermi 
temperature, in which case the finite-temperature contributions extend beyond 
those given by the FLT. In neutron star mergers, it is likely that in some cases a 
hyper-massive neutron star, or HMNS, is formed: the merged remnant mass ex- 
20 ceeds the cold maximum mass. The metastable support is provided by rotation, 
including differential rotation, and thermal effects. The timescale over which 
collapse to a black hole eventually occurs, potentially observable in gravitational 
wave signatures, will therefore be sensitive to thermal effects |1]. In this con¬ 
tribution, we derive analytical formulas for next-to-leading order temperature 
25 effects in the state variables of interacting nucleons in both the non-relativistic 
and relativistic limits for a variety of nuclear interaction models. 

For non-relativistic models with zero-range interactions, knowledge of the 
Landau effective mass is sufficient to satisfy thermodynamic identities. How¬ 
ever, in the general case of finite-range interactions, momentum derivatives of 
30 the Landau effective mass function up to second order are required. We compare 
results from the analytical expressions to exact numerical calculations for zero- 
and finite-range potential models as well as for relativistic mean-field theoretical 
models. The analytic next-to-leading order expressions lead to an improvement 
of the leading order results of FLT, as demonstrated by the wider ranges of de- 
35 generacy and temperature for which they remain valid. In addition, we derive 
relations in a form that are independent of the dimensionality of the momen¬ 
tum space under consideration. Therefore, although our discussion focuses on 
examples from dense matter physics, which are three-dimensional systems in 
momentum space, the expressions derived can have a wider application to cer- 
40 tain problems in condensed matter physics in which the momentum space is 
two-dimensional. 
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The paper is organized as follows. In Sec. the formalism to calculate next- 
to-leading order corrections to the results of FLT in D-dimensions is developed. 
Analytical formulas appropriate for 3-dimensions are given in Sec. whereas 
45 Sec. [^contains results for 2-dimensions. The formalism is applied to zero- and 
finite-range potential models and a relativistic field-theoretical model in Sec. 
Numerical results for these models are presented in Sec. where the extent to 
which the next-to-leading order corrections improve the FLT results are demon¬ 
strated. Section presents a summary and conclusions. Useful formulas for the 
50 evaluation of the thermal properties are provided in Appendices A, B, and C. 


2. GENERAL CONSIDERATIONS 


For a generic Hamiltonian density r) where n and r are the number and 
kinetic energy densities respectively, the single-particle potential U is obtained 
from a functional differentiation of % with respect to n, and can contain terms 
55 depending on n as well as the momentum p\ 

U{n,p)=°^=U{n) + R{p), (1) 


where U{n) denotes contributions that depend on n only. Note that R above 
may also be n-dependent but we will suppress this for notational simplicity. 

The study of the thermodynamic properties of a fermion system involves 
integrals of the form 


I = 


dp g{p)fip), f{p) = 


■ exp 


e-^ 


( 2 ) 


60 where T is the temperature, p is the chemical potential, and the single-particle 
spectrum 

e= ^ + (3) 

The structure of the function g(p) is determined by the state variable under 
consideration. In general, integrals involving the Fermi function f{p) do not 
admit analytical solutions and thus require numerical treatment. In the low- 
65 temperature limit, however, when the degeneracy parameter 


V = 


p-f-{p = 0 ) 

T " 


(4) 


is large, these integrals can be approximately evaluated employing the Sommer- 
feld expansion (see, e.g., 0) by transforming Eq. ([^ to 


1 = 


dy 


Hy) 


1 -h exp(?/ - g) 




d^(j) 

360 dy3 

(5) 
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(6) 


with the identification (j){y)dy = g(j))dp, and the substitution 


y = 

from which it follows that 
dy _ p 
dp 


— U{n) 
T 


2mT 


+ 


R{p) 


MT 


and (^(y) = ^'^g{p ), 


70 where the Landau effective mass function 


M{p) = m 


1 + 


m dR(p) 
p dp 


(7) 


( 8 ) 


This function is implicitly temperature-dependent and its relation to the Landau 
effective mass m* is 

M{p = pf\T = G) = m\ (9) 

where pp is the Fermi momentum. From the relations in Eq. 0 

T^M^g 


d(j) 
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1-p 


M' 

J4 


( 10 ) 
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+ -eP 
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( 11 ) 


75 where the primes denote differentiation with respect to p. 

For a system in D dimensions having 7 internal degrees of freedom, the 
number density is given by 


n = Cd dp p fp with Cd = 


7 


i77r^/2 


(27rV^ {D/2)\ ■ 

The combination of Eqs. ([^,(10), and (11) with = p^~^ results in 


( 12 ) 
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TT 
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Dp^-^MlT^D-2+pp 


Ma 


7^n 

360^^' 

( P,M'p 
\ Mp 


P-S 
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( 13 ) 


+ (6L>2 _ 4011 + 59) 




vlM'M 


Ml 


+ (411-11) 


pIM'I' 
M, 

pIm: 


M, 


where the subscript /i denotes quantities evaluated at e = /x, i.e., 


pI 


e = M = ^ + RiP^i) +U{n). 


(14) 


For N particles in a volume V, the number density n = N/V at T = 0 and at 
finite T is the same. Equating the result in Eq. (131 to its T = 0 counterpart 
n = CdP^p/D, and perturbatively inverting we get 


with 


Pfj. =Pf 


PF = 


1 - 


7r2 


6 


P% 


D-2 


PfM'p 





, . dM 


and Mp = —— 
dp 


(15) 


(16) 


PF 


As our main goal here is to derive the next-to-leading order correction in tem¬ 
perature for the entropy density s, it suffices to truncate the series expansion of 
Pfj, to 0{T^). We will show below that higher-order terms do not contribute at 
this level of approximation where we may also neglect the temperature depen¬ 


dence of Ai and its derivatives. The result in Eq. (151 helps us to work only 


with quantities defined on the Fermi surface as done in Landau’s Fermi-Liquid 
theory [DEIE!. The entropy density is formally given by 


s = -Cd I dp p^ ^ [fp In /p - (1 - /p) ln(l - /p)]. 


(17) 


Integrating this expression twice by parts we obtain 

^ ^ f {m ^ ^ ~ ^ J 

where 


R{p) 


P dR{p) 
D dp 


T = Cd I dp 


(18) 

(19) 


Eq. (18) can be written as 


is the kinetic energy density. With the aid of Eq. (|14|) for the chemical potential, 

P dR{py 



/ dpp^-^fp 


Rip) - RiPfi) 


D dp 


( 20 ) 
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95 from which we identify the functions 


9is{p) = 


P 


92s{p) = P 


m 

D-l 


1 

D 

Rip) - RiPfi) 


-1 Pu 


^D-1 


2m 

P dRjp) 
D dp 


( 21 ) 

( 22 ) 


to be used in the Sommerfeld expansion. For both of these functions, the hrst 
term on the right-hand side of Eq. ([^ involving an integral vanishes yielding 


^ +^-ZrCDP^-''MlT^ 


-Copi 


4/2 


.M': 


90 


M ■' // 

M', 


(D-4)(D-2) 




(23) 


Use of Eqs. (15) and (16) in the above result delivers the working expression 


for s in terms of quantities dehned on the Fermi surface: 


s ~ 




(D-9)(D-2) 


7 

2 


7 


2 m* 2 TO* 


(24) 


where the 0{T) term is the well known result from FLT. We note that a large 
number of cancellations occur in obtaining Eqs. (23) and (p4|) despite the com¬ 


plexity of of Eqs. (10) and (11). For a system composed of different kinds of 


particles the total entropy density is a sum of the contributions from the indi¬ 


vidual species where, in Eq. (24), the Fermi momentum, the effective mass, and 


105 its derivatives all carry a particle-species index i. 


Equation (24) forms the basis from which other properties of the system can 


be derived. Eor example, the entropy per particle is the simple ratio S = s/n, 
whereas the thermal energy, pressure, and chemical potential can be obtained 
through the application of the appropriate Maxwell relations [S] : 


Eth= TdS, Pth = -n^ 




and 


Pth = 


an 


(25) 


110 [for a multiple-species system, pth,/ = — Jids/dni)dT]. 

The specihc heats at constant volume and pressure are given by the standard 
thermodynamics expressions [5] 


Cp 


= T 


= T 


95 

dn 

95 

dn 


dEth 

dT 


= Cv 


rji ^ 9Pth 


dT 


dP l 
dn It 


(26) 

(27) 
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We note that the formalism above has not considered effects, for example, 
from single particle-hole excitations, or from collective and pairing correlations 
115 near the Fermi surface Contributions from these sources must be added 

to those considered here when appropriate. 


3. RESULTS FOR D=3 


For a single-species system of spin 1/2 particles in 3 dimensions [for which 
C 3 = l/(7r^?i^)], the entropy density becomes 


s = 


PFm*T 27r2 m*^T^ .. ^ , 




120 where 


Lf = 


M 


12 ^^ 771*2 


'2 7 2 3 

12 777 * 4 


M' 


We stress that, in general, 

dM{p) 


dp 


= M'f ^ TO*' = 


dM{pF) 

dpF 


(28) 


(29) 


(30) 


as R can contain both p and pf (via n). In terms of the level-density parameter 
a = 7r^7T7*/(2p|.) = 7r^/(4TF) (where Tp is the Fermi temperature), Eq. (28l 
can be written as 


16 

s = 2anT - -a^nT^il — Lf)- 

0772 


(31) 


125 The quantity Lf arises from nontrivial momentum dependencies in the single¬ 
particle potential. For free gases (where R{p) = 0), and for systems having only 
contact interactions where R{p) oc p^ (such as Skyrme models), Lf = 0- 
Equation (31) in conjunction with Eqs. (25 )-([^ leads to 

S' 


2aT - - Lf) , Eth = aT^ - - Lf) (32) 


Eth 

Pth 

f^th 

Cv 

where 


aT'^ - - Lf) 

OTT^ 

‘^anQT‘^ - -^a^nQT^ 
3 OTT^ 


1 — Lf 


n dLF\ 


2Q dn 




2Q 




(1 - Lf))! - 2Q) - n 


dL p 


dn 

16 


40 

and Cp = Cv + g ^ 


dn 


Q = l- 


3n dm* 
2 m* dn 


(33) 

(34) 

(35) 

(36) 

(37) 


130 In the derivation of Eq. (36) we have assumed that the zero-temperature pres¬ 
sure Po is such that dPo/dn 3> dPth/dn. This condition will not be met in 
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situations where Pq is relatively flat as in the vicinity of a critical point. When 
this is the case, we must use Eq. (271 for Cp, with 


fdPth 


V dT 


^anQTj 


dP 

dn 


4 

3 

dPo 

dn 


1 48 2 rp 2 (^ r I ^ dLp 

57 r2“ (, 2Q dn 

2Q dQ 


(38) 


+ —aQT 1 — ——h 


dn 




l-2Q + n^ ) ( 1-Ef 


n dLp 
2Q dn 


Tl —;- I 1 — 2Q + 


dn 


n^ d^Lp 

2 Q'^ dn J ' 2 Q dn^ 


+ 


(39) 


Similar considerations as with Eq. (pM|) hold for the ratio of the specihc heats 

(40) 


Cp 8 aQ^T^ 


dn 


135 Other quantities of interest in astrophysical applications include the thermal 
index Tth 


Lth — 1 + 
and the adiabatic index E 5 


Pth 

nEth 


= 1 + ?^- 


57r2 


a^nT 


, dLp 
dn 


_ C^n ^ 


n 

'K 


dPo 


dn 


2Q dQ 


dPn 


^ + -aQT^ 1 + ^ + 


dn Pq dn 


where, in addition to Eqs. (34),(39l and (40), the approximation 


1 


1 




_ _ Pth 

p~yo ''" 

was used. In its native variables (n, S), Fs is given by 

n 


n dP 
P 


Pn 


nQS^ 

6a 


dPp 

dn 


6 a V 3^ 


n dQ 
Q dn 


(41) 


(42) 


(43) 


(44) 


140 To arrive to Eq. (|44|) one begins by inverting Eq. (p^ for the small parameter 

(45) 




which is then employed in the expression for the thermal pressure with the 
results 


Pu. = 

oa 


nQ ^4 


dOTT^a 


1 — Lp — 


3n dLp 
2Q dn 


(46) 


dPth 


dn 


Q 


6 a 


n dQ 


= 1+oQ+7^^ 


Q dn 






































(47) 


Q 


— 2n 


SOTT^a 
dL p 


5" 


dn 


2 n dQ\ 

' + 3« +<3 

3 \ 3n^ d^Lp 

~ W dn^ 


Finally, the result is truncated to 0{S^) in both the numerator as well as the de¬ 
nominator. We refrain from invoking approximation (43) as for nuclear systems, 


Pq can cross 0 at low densities. This is not a problem in the variables (n, T) 
because the degenerate approximation breaks down at sufficiently low density 
regardless of T. In the variables (n, S), however, for small values of the entropy 
the system remains degenerate irrespective of the density, and thus division by 


zero is avoided (as could happen if Eq. (43) is used). 


We point out that the adiabatic index is related to the squared speed of 
sound Cs according to 

= (48) 

V c/ a + mn 

where h = nE + P is the enthalpy density. 


4. RESULTS FOR D=2 


In condensed matter physics, 2-dimensional systems are of much interest. In 
155 the current framework, the entropy density is 


2 






s--C2m T+—C2^-ypp^+Pp^-pp-^^ 


(49) 


with C 2 = (l/27r7i^). A noteworthy feature of this result is that the 0(T^) term 
receives contributions only from the derivatives of the effective mass function 
with respect to p at the Fermi surface. Thus, it is absent not only for free gases 
but also for systems with contact interactions where the p^-dependence of R 
160 implies that dM/dp = 0. 

In terms of the level density parameter a = 7r^TO*/(2p|n), and 

1/2 

, Q = l- 


Pf = 


2 n\ 


dm* 


Lp = 


Pf 


M', 


Eq. (49) leads to 


S' 

Pth 

Pth 

Cv 


= laT- 


56 

457r^ 


m 


a^T^Lp , 


-Pf- 


TO’ 


-PF- 


m* dn 

M'p 




= \anQT‘^ + 

3 lOTT^ 




Lp — 


n dLp 
3Q dn 


= 


56 

I57r^ 


a^T^Lp 


and 


Lp{l-‘^Q)+n 

Cp = Cv + 


dL 1 


dn 
16 

q dPp 

dn 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 
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The above results do not include the effects of collective excitations near the 
Fermi surface or of non-analytic contributions. As pointed out in Ref. 

165 2-dimensional Fermi systems in condensed matter physics (even with contact 
interactions) have contributions to the entropy from interactions separate 
from those due to the collective modes. These contributions arise from non- 
analytic corrections to the real part of the self-energy. 


5. APPLICATION TO MODELS 

170 In what follows, we compare the analytical results from the leading order 
corrections to Landau Fermi-liquid theory to the results of exact numerical cal¬ 
culations of the thermal state variables. These comparisons are made using 
models that are widely used in nuclear and neutron star phenomenology. In the 
category of non-relativistic potential models, we begin with the model, referred 
175 to as MDI(A), that reproduces the empirical properties of isospin symmetric 
and asymmetric bulk nuclear matter cni, optical model fits to nucleon-nucleus 
scattering data m, heavy-ion flow data in the energy range 0.5-2 GeV/A [T2] . 
and the largest well-measured neutron star mass of 2 Mq [Hill]. This model, 
which is based on Refs. [HKle], incorporates finite range interactions through 
180 a Yukawa-type, finite-range force, is contrasted with a conventional zero-range 
Skyrme model known as SkO' m- Both models predict nearly identical zero- 
temperature properties at all densities and proton fractions, including the neu¬ 
tron star maximum mass, but differ in their predictions for heavy-ion flow data 
[Is] . To provide a contrast, we also investigate a relativistic mean-field theo- 
185 retical (MFT) model [T0| which yields zero-temperature properties similar to 
those of the two non-relativistic models chosen here. For all three models, we 
consider nucleonic matter in its pure neutron-matter (PNM, with 7 = 2 ) and 
symmetric nuclear matter (SNM, with 7 = 4 ) configurations. 


5.1. Finite-range potential models 

190 For the MDI(A) model [111 HD] j the momentum-dependent part of the single¬ 
particle potential is given by 


R{p) = 


2 C^ 


d^p' /, 


no (27r?i)3 


p' 


T=o a A3 


no TT^hi^ 
(A^+pI-p^) 
4Ap 


^ — arctan 

A 


i + (*¥) 

P+Pf 


A 


In 


A^ + ip + ppf 


(56) 


arctan 


P-PF 

A 


_A 2 {p- pfY_ 


(57) 


For the coefficients Uq, C 2 , C '4 and A we use the values 0.16 fm“3^ -23.06 MeV, - 
128.9 MeV and 420.9 MeV, respectively. Explicit expressions for the derivatives 
of R{p) and their connection with Ai and Lp are provided in Appendix A. 
195 The MDI Hamiltonian density is shown in Appendix B. For details of the exact 
numerical calculations, see Ref. [TO] . 
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5.2. Zero-range Skyrme models 

Zero-range Skyrme models belong to that subset of the D — 3 case for which 
Lp = 0. This is because, for these models, the momentum-dependent part of 
200 the potential has the form 

R = I3{n)p^ (58) 

{P{n) is a density dependent factor) which renders the generalized effective mass 
to be p—independent: 

777 1 

“ 1 + “ 1 -h 2mf3{n) ’ 

and therefore its derivatives Ai' = M." = 0. Hence Lf = 0 as well. Conse¬ 
quently, the results in Sec. [^for Skyrme models simplify considerably. Results 
205 to be shown here are for the SKO' model ini, the exact numerical calculations 
for which are described in Ref. m- 


5.3. Relativistic models 

The single-particle energy spectrum of relativistic mean-field theoretical 
models m obtained from the nucleon equation of motion has the structure 

e = E* +U{n), E* = [p^ + . (60) 


The single-particle potential U (n) is the result of vector meson exchanges whereas 
the Dirac effective mass M* arises from scalar meson interactions. The imple¬ 
mentation of the above equations in the Sommerfeld expansion is made possible 
by the identification 


y = 


E* 

Y 


Y 

dp 


P 


E*T 


E*T 

and (j){y) = - g{p). 


(61) 


The calculation of d(j)/dy, d^(j)/dy^ and n proceeds as in the non-relativistic case 
with the replacement M. ^ E* [cf. Eq. ®]. In particular, for p^ we have 


Pn =Pf 


1 - 


EfY 

P% 


D-2 


PpEf 

E% 


(62) 


where 


E*p = Yf + 


*2 


)i/2 and 

dp 


PF 

E* 


(63) 


The simple dependence of E* on the momentum p in Eq. (60) leads to the 


correspondingly straightforward expression (63) for Ep which, as we will show 


soon hereafter, results in an elementary form for Lp and by extension the whole 
set of the MET thermodynamics can be written in an uncomplicated manner. 


Substituting Eq. (63) into Eq. (62) yields 


Pn =PF 


1 - 


TT^ E*JY 


6 


P% 


D-2 


Pf 

Ef 


(64) 
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The twice-by-parts integration of Eq. 0 for the entropy density in the rela¬ 
tivistic context gives 


C. 


s=^ I dpp^-^U{E* + 


P dE* _ 

D dp ^ 


(65) 


where one observes the analogy with the integral term of Eq. (20). Using 

p dE* 


9s{p) =P 


D-1 




77'* 

D dp 


( 66 ) 


225 we proceed as before to get the entropy density in terms of p^ as 


Ttt 


= ^CdP^-'^E*T + ^Cd{D - 2){D - A)p^-<^E:^T^ 


1 + 


{D - 4) Ef 


(67) 


which, with the aid of Eq. (62), becomes 
^2 


S = '—CDPF~'"E*pT+'—CD{D-2){D-'d)p^-^E*p^T^ 


11 


pI 


Pf 


( 68 ) 


2{D-9)E*p^ 2{D - 2){D - 9) E*p^_ 

In the derivation of the last equation the weak temperature of M* in the de¬ 
generate limit has been ignored (but not of E*). Combining Eq. ( 68 ) with 
Eqs. (|2^-(27) in H = 3, and using the definitions [here, the Fermi temperature 

230 Tp = Pf/{2E*p)] 


TT^ E*p 
2 Pf 

we obtain 


q = 


M *2 


1 - 


3n dM* 
M* dn 


j r Pf ^ Pf 

and Lf = - 


12 Et 


(69) 


5 = 2aT- ^a^T^l-Lp), Eth = aT^ - - Lp) (70) 


Pth = ^anT^(l + q)- ■^a^nT‘^ 


T _ r I ( 1 _ hn _ 


3 9 E*p^ J 




(71) 

(72) 


Cv = 2aT- -^a^T^il-Lp) and Cp = Cy + ^. ( 73 ) 
OTT^ 9 


dn 


As in the nonrelativistic case, when conditions are such that dPo/dn is small, 
one must use derivatives of the pressure with respect to n and T that include 
thermal contributions to 0{T'^) in the calculation of Cp. Explicitly, 


(dPth 
^ dT 


= -a^n^T‘^{l + qf 
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dP 

dn 


T 


dq 

dn 


1 - 

dPo 

dn 

4 

5tt^ 


48 


Stt^ 1 
T^a 


■ q 


1 — Lp P qil — 




(1 + q){2 - q) + 3n 


dq 

dn 


1 — Lp 


rr^A 3 

1 a q 

n dLp 
q dn 

2 , ,1 dM* 

-<z) + 


n dq 
q dn 


1 - 


10 ^ 

9 Ep 


fl+9\ 40^ 

V ^ sM 27 E*P 


( 75 ) 


M* fdM* 3nd^M*\ 


235 6. RESULTS AND DISCUSSION 

Here, we compare the results from FLT and FLT+NLO with the exact nu¬ 
merical results for the two non-relativistic models (MDI(A) and SkO'), and for 
a relativistic mean-held theoretical model (MFT). Numerical techniques for ob¬ 
taining the exact numerical results are detailed in Refs. nnun]. The thermal 
240 properties presented for PNM and SNM are at a temperature of T = 20 MeV. 


In the top panels of Fig. the Landau effective masses m* of the neutron 
scaled with its vacuum value are shown as a function of baryon density n. For 
the MFT model, both m* = Ep = ^M*^ + and the Dirac effective mass 

245 M* are shown. Noteworthy points for the non-relativistic models are: (i) The 
isospin splittings are qualitatively similar - m*/m being larger for PNM than 
for SNM - although quantitative differences are present, and (ii) except for n 
up to 0.2 fm“^, the decrease with increasing n for the MDI model is relatively 
slow (logarithmic decline) compared with that for the SkO' model [(1 -I- /3n)“^ 
250 fall off]. This overall flatness of m* for the MDI model is a direct consequence 
of the momentum structure of its single-particle potential which causes it to 
saturate at high momenta. For the MFT model, M* decreases monotonically 
with n to values lower than those for the non-relativistic models. The Landau 
mass m*, however, exhibits a non-monotonic behavior, attaining a minimum for 
255 nmin = 0.57 (0.52) fm~^ obtained from the solution of 


Pp dM* 
M* dpp 


(77) 


for the case of SNM (PNM), and increasing thereafter due to the monotonic 
increase of kp^ with n. Physically, nmin marks the transition of nucleons well 
into the relativistic region. The density np at which pp = M* occurs at 

/M*\^ _o 

np = 0.643 7 f -j fm (78) 
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about (2/3) Umin = 0.38 (0.34) fm ^ for SNM (PNM) which signals the onset 
260 of relativistic effects which become progressively important for n > ur- 

The bottom panels of Fig. show the logarithmic derivatives of m* vs n. 
Also shown for the MFT model is dXnM*/dhin which has been divided by a 
factor of 3 to fit within the figure. The logarithmic derivatives m* for MDI(A) 
265 show little variation with n at supra-nuclear densities. In contrast, results for the 
SkO' model, which take the simple form (m*/m) —1, show a significant variation 
with n. The logarithmic derivative of M* in MFT drops to values considerably 
lower than for the other two models. This derivative remains negative but ap¬ 
proaches a constant value at large densities. The logarithmic derivative of the 
270 Landau effective mass drops until ^ 2no) but then increases to positive values. 
This behavior is a reflection of the minimum that occurs for m* in this model. 
As will be seen below, the density dependences of the effective masses and their 
logarithmic derivatives determine the behavior of all the thermal properties in 
FLT. Higher order derivatives of the Landau effective mass function in Eq. 

275 appear in FLT-I-NLO. 

The FLT and FLT-I-NLO results for the thermal energy, Eth vs n, are com¬ 
pared with the exact numerical results in Fig. The NLO corrections to FLT 
yield agreement with the exact results down to sub-nuclear densities of 0.5 to 1 
280 no compared to 2-3 n-o for FLT. As is the case with FLT, slightly but system¬ 
atically better agreement with FLT-l-NLO occurs for PNM than for SNM for 
all the thermodynamic quantities we study. This is a consequence of the fact 
that the neutron density in PNM (n„ = n) is twice the neutron density in SNM 
(n„ = n/2); PNM is more degenerate than SNM at the same baryon density n. 

285 

In Figure we show the convergence of FLT and FLT-I-NLO results to 
the exact numerical results for Pth- For all three models, the FLT-I-NLO re¬ 
sults extend the agreement with the exact results to lower densities than those 
of FLT. The SkO' model shows the greatest changes relative to FLT for both 
290 SNM and PNM. Both the MDI(A) and MFT results also show improvement 
using FLT-I-NLO with better agreement occurring for PNM. For Pth at n > no, 
the influence of m*(n) and its logarithmic derivative with respect to n (in FLT), 
and the higher derivatives of AI (in FLT-I-NLO) are amply demonstrated: (i) 
for MDI(A), Pth grows very slowly with n; (ii) considerably larger growth with n 
295 is exhibited for SkO' than for MDI(A); and (iii) for MFT, the prominent peaks 
in Pth at intermediate densities are due to the minima in m* and din M* /dlnn 
at similar densities. For asymptotic densities, Pth oc n^/^ in MFT models char¬ 
acteristic of massless particles. A corresponding behavior is also present in the 
thermal chemical potentials (see below). 

300 

Comparisons of the entropy per baryon, S, are shown in Fig. For all 
three models, the agreement between the exact results and those of FLT-I-NLO 
extends to below noj and to as low as no/4 for PNM with MFT. This is an 
improvement from the FLT results for which convergence ranged from 1.5no 
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(FLT and PNM) to 3no (SkO' and SNM). 


In Fig. the thermal parts of the neutron chemical potentials, Hn,thi are 
shown as a function of n. Results from FLT+NLO lie closer to the exact results 
than do those of FLT with the agreement extending to sub-nuclear densities. As 
310 with the other state variables, the agreement is quantitatively better for PNM 
than for SNM. 


The specific heat at constant volume, Cy, is presented in Fig. |^as a func¬ 
tion of n. For all models, the FLT-I-NLO results show better agreement with 
315 the exact results than those of FLT. The lowest density for which the agreement 
extends differs between the models, that for the SkO' model being higher than 
for the other two models particularly for SNM. 


In Figure we show the specific heat at constant pressure and its limiting 
320 cases as functions of n. The maxima at low sub-nuclear densities in the ex¬ 
act numerical results are related to the liquid-gas phase transition of nucleonic 
matter which occurs at T ~ 15-20 MeV. The convergence between the exact 
and approximate results for MDI(A) and MFT follows a pattern similar to the 
quantities discussed previously with the conspicuous exception of FLT which 
325 appears to outperform FLT-I-NLO for SNM. We attribute this feature to a nu¬ 
merical accident, possibly due to the proximity to the phase transition. Note 
that FLT begins to deviate from the exact result for densities below n ^ 3no in 
the case of SkO'. As our analysis here is concerned with the degenerate region, 
this failure at sub-nuclear densities is not surprising. For an adequate treatment 
330 in the non-degenerate region, see Ref. m- 


The thermal index, Fj/j = 1 -|- (Pth/^th), is shown in Fig. For the SkO' 
model, the FLT result is exact for all regions of degeneracy. This happenstance is 
due to the fact that for non-relativistic nucleons with only contact interactions, 
335 Pfii and Eth can be written entirely in terms of their ideal-gas counterparts as 


Pth{n,T) = Pl^{n,T-m*)Q 

-pid 9 

eth{n.T) = - (79) 


340 


regardless of the degree of degeneracy. The n-dependence of m* for Skyrme-like 
interactions thus yields 


Pt?i 


8 m*{n) 

3 m 


(80) 


For the finite-range MDI(A) and the relativistic MFT models, the FLT-I-NLO 
results mildly improve the FLT results in reproducing the exact ones. This 
marginal improvement is related to the ratio Pth/sth in ^th in these cases, 
which amount to a constant plus a correction due to effects of n plus a second 
correction due to (n, T) effects, which means that temperature effects are sub¬ 
leading, and therefore very weak. Note that the results for the three models 
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differ significantly from each other, both qualitatively and quantitatively. These 
345 differences are due to the differences in the effective mass functions of the models. 

In physical applications involving neutron stars and supernovae, contribu¬ 
tions from leptons (electrons and muons) and photons must be included to all of 
the state variables. Thermal effects from these sources are adequately given by 
their free gas forms, and numerical methods for their calculation for arbitrary 
350 degeneracy can be found in Ref. [20]. The influence of thermal effects from 
leptons and photons to the total has been detailed in Ref. [10], and will not be 
repeated here. 

7. SUMMARY AND CONCLUSIONS 

For homogeneous systems of fermions in the limit of extreme degeneracy 
355 (T/Tp ^ 1, where Tp is the Fermi temperature), Landau’s Fermi Liquid The¬ 
ory (FLT) provides simple analytical expressions that are model independent 
for the thermal state variables {e.g., entropy, energy, pressure, chemical poten¬ 
tial, and specific heats) [T]. In the absence of collective excitations close to the 
Fermi surface, thermal effects are primarily determined by the nucleon’s Lan- 
360 dau effective mass and its hrst density derivative which in turn depend on the 
momentum-dependence of the T — 0 single-particle energy spectrum. 

In this work, we have developed a method by which thermal effects in near¬ 
degenerate to degenerate matter can be described to next-to-leading order in 
T/Tp for models with general momentum dependences in their single-particle 
365 potentials. Analytical formulas valid to next-to-leading order in T/Tp for all 
of the thermal state variables are presented. The entropy density and spe¬ 
cific heats are carried to 0{T/Tp)^ whereas the energy density and pressure 
to 0{T/Tp)‘^, extending the leading order results of FLT. These extensions in¬ 
volved the use of a generalized Landau effective mass function which enables 
370 the calculation of the entropy density, and thereafter the other state variables, 
for a general single-particle spectrum. In special cases, e.g., models with con¬ 
tact interactions, knowledge of the Landau effective mass suffices. In the case 
of finite-range interactions, momentum derivatives of the Landau effective mass 
function up to second order are required to satisfy the thermodynamic identity. 
375 Our results are valid for potential and field-theoretical models as long as the un¬ 
derlying interactions yield a single-particle spectrum that is weakly dependent 
on temperature in the degenerate limit. 

We find that 0{T^) corrections in 2-dimensions for non-relativistic models 
appear only if they include finite-range interactions or, equivalently, if the mo- 
380 mentum content of their mean field is something other than quadratic. These 
contributions supplement the O(T^) non-analytic contributions even with zero- 
range interactions established previously [9j. 

To illustrate the density region of their applicability, numerical results from 
the new formulas were compared with those of exact numerical calculations 
385 for zero- and finite-range potential models as well as for relativistic mean-field 
theoretical (MFT) models widely used in astrophysical applications of hot and 
dense nuclear matter. In all cases, excellent agreement with the exact results 
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was found even to sub-nuclear densities of ^ 0.1 fm“^ for T = 20 MeV, whereas 
FLT results are valid only for densities beyond ^ 0.3 fm“^. (For low tempera- 
390 tures and below ~ 0.1 fm~^, inhomogeneous phases with nuclei and pasta-like 
configurations are known to exist, and must be treated separately.) In addition 
to providing physical insights, our analytical results facilitate a rapid evaluation 
of the EOS in the homogeneous phase (important for computer-time consuming 
large-scale simulations of supernovae, proto-neutron stars, and mergers of binary 
395 compact objects). Our formulas can be used for any T = 0 quasi-particle spec¬ 
trum, e.g., those extracted from Brueckner-Hartree-Fock and Dirac-Brueckner- 
Hartree-Fock approaches, extensions of MFT models with non-linear derivatives 
or with 2 -loop effects, and effective field-theoretical approaches. 

Examples of contributions not included in our work arise from, e.g., non- 
400 analytic contributions from single particle-hole excitations and, collective and 
paring correlations close to the Fermi surface. Their roles as functions of den¬ 
sities and temperatures of relevance to astrophysical phenomena need further 
investigation. 
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Figure 1: Top panels: Effective mass vs density for non-relativistic potential models (MDI(a) 
and SkO^) and relativistic mean-field theoretical model (MFT) for symmetric nuclear mat¬ 
ter (SNM) and pure neutron matter (PNM). Bottom panels: Logarithmic derivative of the 
effective mass with respect to density vs density. 


18 










0 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

n 


Figure 2: Thermal energy vs baryon density for the three models at a temperature of T = 20 
MeV. Results for SNM are in the top panels and for PNM in the bottom panels. 
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Figure 3: Same as Fig. but for thermal pressure. 
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Figure 4: Same as Fig. but for entropy per baryon. 



Figure 5: Same as Fig. but for thermal neutron chemical potential. 
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Figure 6: Same as Fig. but for specific heat at constant volume. 
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Figure 7: Same as Fig. but for specific heat at constant pressure. 
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Figure 8: Same as Fig. but for the thermal index. 
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Appendix A. RESULTS RELATED TO Lp 

Here we collect relations required to calculate the quantity Lp defined in 
Eq. @ of Sec. § 


410 Model-independent relations 


M' = 


M" = 


dLp 

dn 

dLp 

dpF 


dM 

dp 

2 M^ 


fdR 


d^R 


dp ^ dp^ 


P 

Pp dLp 
3n dpF 

5 M'p 
12 m* 
m* d^R 

6 dp^ 


pz 

dR 

dp ^ dp"^ 

1 -^ 1 - 
P 


d‘^R\ 2M^ fdR d'^R , p^ d^R 

2 dp^ 


pO 

p dp^ dp'^ 


dp ^ dp'^ 
d'^M 


4 M'p 

^-PF -- 

m* 


, 42 M'p 

1 - -rPF —r 
5 m* 


Pf 


49 M'p 

1+ -R-PF —T 
2 m* 


7 


d d^R 


12 ^^"^ dpF dp^ 


The subscript F denotes evaluation at p = pf- 
MDI model-specific realations 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 


Pf 


dR 

dp 

d'^R 

dp'^ 


d?R 

dp^ 


PF 


d d^R 

dp p dp^ 


A 2 


no TT 




'- 2 '' 


2 C^ M 


2 , In ( 1 + ^ 
2 plJ \ M 


(3p|+A2) 


A _2 


no ppiipj, + A'^) . 

{p%+A^) (1 + ^) 

2 C^ 1 1 

no p|(4p|, + A2)2 

8 p% + 34p|,A^ + 21p|A^ + 3A® 


~r ( 1 + IT ) (‘^Pf + In ( 1 + 


3A2 


Pi 


A 2 


2(7^ 1 


3A^{p% + A^){ip% + A^)^ In ( 1 + ^ 


(A. 6 ) 


(A.7) 


(A. 8 ) 


no p'^h^ Ppiipl + A2)3 

4p%{24p% + 134p|A2 + llSploA^ + 33p|A® + 3A®) 


(A.9) 
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Appendix B. MDI HAMILTONIAN DENSITY 


For the MDI models, the Hamiltonian density is composed of terms arising 
415 from kinetic sources, Hk, density-dependent interactions, Hd, and momentum- 
dependent interactions, 'Hm'- 


% = %k + Hd + 'Hm ■ 


(B.l) 


At T = 0, 
Hk 
Hd 
-Hm 

with 


1 , , _ 1 1 


(PFn +PFp) 


(B.2) 


+ ^n^{l - 2xf + [1 - 2/(1 - 2xf] (B.3) 

2no 2no cr + 1 n, ■' 


0 


^(i +/ ) + ^/ 

_ \^nn \ -^pp) ' ^np 


no 


X 


J-ij 


no 


Up/n, ppi = 

8 ^ 2 a 2 r , 

A A 


2 X PFiPFj^? 
PFj) ^ 


4A, 3 
-Y{PF^■ 

4A, 3 
-J{PF^■ 


Pp^) arctan 

-p?.,.)arctan(^P^i±^' 


A 


A4 , {p%+pj,^)A^ 

-- - - 

12 2 


1-^2 


(B.4) 


(B.5) 


(B. 6 ) 


In this work we use the coefficients Ai = —69.48 MeV, A 2 = —29.22 MeV, 
420 B = 100.1 MeV, cr = 1.362, y = -0.0328, Ci = -23.06 MeV, C7„ = -105.9 
MeV, and A = 420.9 MeV [TU]. 

When a calculation for multiple-species is undertaken, highly asymmetric 
configurations should be avoided as the various particle types involved will be in 
different regimes of degeneracy thus resulting in a slower convergence relative to 
425 the single-species case. Furthermore, one should also refrain from using multiple 
species results with cc = 0 or 1 as numerical complications arise; namely, division 
by 0 occurs in the level density parameters a^. 


Appendix C. MEAN-FIELD THEORETICAL MODEL (MFT) 


The MFT model used here involves the exchange of a, uj and p mesons 
430 (scalar, vector and iso-vector, respectively) [19]. Its Lagrangian density is 


C 




- ToffcoWo - 7oyPoFs - (M 


9<y^o) 
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which yields the following meson equations of motion: 



= g^ns = mlao + '^glal + 

(C.2) 

Wo 

= ^(d/tvl/) = 

(C.3) 


mi mi 


Po 


(C.4) 


in the mean-field approximation with classical expectation values denoted by 
the subscript “0”. The equation of motion for the nucleon ('!') held is 

-f yPoT-s) - M* = 0 (C.5) 

where M* = M — and yields the nucleon single-particle energy spectrum 


e,± = ±El + - n ,), E* = {p} + . 


435 The subscripts z, j refer to the nucleon species, the positive sign to the particles 
and the negative sign to the antiparticles. The thermodynamics of the system 
are obtained from its energy-momentum tensor 


= + ^ rnlal ++ ^{ga(ro)^ - mlwl - mlpl . 

(C.8) 


For an isotropic system in its rest-frame, the energy density and the pressure 
are given by the diagonal elements of T^i, as 

+ ~+ b(ffcrO’o)^ + — (ffcrCro)"' (C.9) 


SrrZp ^ 


p = = f f, 


■ pj I 9l 2 

(pI + (27r?z)3 2m2, 




+ ^ rnlal +-{g^aof +—ig^ao)"^ ■ (C.IO) 


440 The minimization of the grand potential D, = —PV with respect to cto (equiva¬ 
lent to de/dao = 0 at T = 0 and to dP/dao = 0 at hnite temperature) leads to 
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(C.ll) 


a self-consistent equation for the Dirac effective mass 


M* = M - 


mf, 


Us - ^{M - M*f -^{M - M*f 
2 6 


In the present work we use the masses M = 939.0 MeV, ma = 511.2 MeV, 
= 783.0 MeV, rup = 770.0 MeV and the couplings = 9.061, = 10.55, 

445 gp = 7.475, At = 9.194 MeV, A = —3.280 x 10“^. These correspond to a cold 
symmetric nuclear matter equilibrium density no = 0.155 fm“^ at which the 
energy per particle E/A = —IQ MeV, the compression modulus Kq = 225 MeV, 
and the symmetry energy Sy = 30 MeV. 
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